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The rheological properties of soft materials often exhibit surprisingly universal linear and non-linear 
features. Here we show that these properties can be unified by considering the effect of the strain- 
rate amplitude on the structural relaxation of the material. We present a new form of oscillatory 
rheology, Strain-Rate Frequency Superposition (SRFS), where the strain-rate amplitude is fixed as 
the frequency is varied. We show that SRFS can isolate the response due to structural relaxation, 
even when it occurs at frequencies too low to be accessible with standard techniques. 

PACS numbers: 83.60.Bc, 83.60.Of, 83.80.Hj, 83.80.Iz 



Concentrated suspensions [H , pastes 0], emulsions [1], 
foams [1] , or associative polymer systems [1, 0] represent 
commonly encountered examples of soft solids. Their 
structures are often metastable, with dynamics strongly 
reminescent of that of glasses. They are thus often 
termed soft glassy materials Q . Despite the wide di- 
versity of these materials, their rheological response, as 
characterized by the complex shear modulus G*(u>) — 
G' (uj) + iG" (uj) , is remarkably similar. The storage mod- 
ulus G' is only very weakly dependent on frequency. The 
loss modulus G" also shows a very weak frequency de- 
pendence and is much larger than the response due to 
the continuous phase fluid. It often exhibits a shallow 
minimum near the lowest experimentally accessible fre- 
quencies. At even lower frequencies, there should be a 
crossover from solid-like to liquid-like behavior which is 
signaled by a pronounced peak in G"(ui), reflecting the 
structural relaxation. Unfortunately, the relaxation fre- 
quencies are often much too low to be accessed by stan- 
dard linear rheological measurements. A typical example 
of this linear viscoelastic behavior is shown in Fig. 1(a) 
for soft hydrogel spheres. Remarkably, the similarities 
in the rheological response of these materials extend 
even to nonlinear measurements, characterized by strain- 
dependent viscoelastic measurements performed at con- 
stant uj, varying the strain amplitude. Above a critical 
yield strain, the storage modulus exhibits a power-law 
decay G'(-fo) oc Jq U . By contrast, the loss modulus ex- 
hibits a well-defined peak, before falling as a power-law 
G"(7o) oc typically with v" w v'/2. The pro- 

nounced peak in the loss modulus is a remarkably robust 
feature of soft glassy materials [E S IE 0] • The ubiqui- 
tousness and similarity of the rheological response, both 
linear and nonlinear, of so many soft materials suggests 
that the response is governed by a common underlying 
mechanism. A possible clue to the origin of this behavior 
comes from a proposed explanation for the response of 
a supercooled fluid [t| . Within this picture, the peak in 



G"(jo) observed in nonlinear measurements is directly 
related to a decrease of the structural relaxation time 
with increasing shear rate [1, U, 12, 13]; the ap- 
plied strain drives the relaxation and forces it to a higher 
frequency, where it is directly probed. While this pic- 
ture provides an excellent description of a colloidal su- 
percooled liquid, the underlying physical concept should 
be much more generally applicable. If this link is veri- 
fied, the relaxation could be probed using the nonlinear 
response, even when the relaxation frequency is experi- 
mentally inaccessible. This would provide a new probe 
of the dynamics and rheology of soft materials. However, 
there have been no attempts to exploit this behavior and 
explore its general applicability to a variety of different 
materials. 

In this Letter, we show that the typical linear and 
nonlinear viscoelastic behavior of a wide variety of soft 
materials can be explained by a common mechanism 
d, [ll , ID, 13j • We introduce a new approach to oscil- 
latory rheology; by keeping the strain rate constant dur- 
ing the measurement we can directly probe the structural 
response. The data exhibits a characteristic scaling, al- 
lowing us to clearly isolate and probe the response due to 
structural relaxation. This new technique, which we call 
Strain-Rate Frequency Superposition (SRFS), should be 
generally applicable to the study of many soft materials. 

To elucidate the viscoelastic response, we use a sys- 
tem of hydrogel microspheres dispersed in water. They 
consist of interpenetrating polymer networks of cross- 
linked poly N-isopropylacrylamide and polyacrylic acid 
14 , [l5| . The sample studied has a polymer weight frac- 



tion c p = 0.0186 and at low volume fraction the particles 
have a radius of R ~ 95nm. All rheological measure- 
ments are performed with a strain-controlled rheometer 
(ARES), applying an oscillatory strain 7(f) to the sample 
and measuring the resultant time-dependent stress cr(t). 

This system of soft spheres shows the behavior typ- 
ical of soft materials, both in its linear and nonlinear 
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FIG. 1: Oscillatory measurements of G'(solid circles) and 
G"(open circles) for a suspension of hydrogel particles, 
(a) Linear frequency dependent measurement at 70 = 0.02. 
The dashed line is proportional to y/uj. (b) Strain dependent 
nonlinear measurement at ui = 1 rad/s. 



viscoelastic behavior, including the ubiquitous peak in 
G"(7o) that typefies the rheology. The measurements 
shown in Fig. 1 represent the most commonly used types 
of oscillatory rheological measurements. They depend on 
two variables only, the frequency u> and the strain am- 
plitude 70; the applied oscillatory strain for each data 
point is 7(4) = 7osin(w£). Frequency dependent mea- 
surements [Fig. 1(a)] are usually performed in the linear 
viscoelastic regime at a constant 70 <C 1, while strain de- 
pendent measurements [Fig. 1(b)] are perfomed at con- 
stant u>, probing the viscoelastic response at a fixed time 
scale 2it/uj. However, the probe frequency might not be 
the only important time scale. Instead, as the strain rate 
becomes large, it can itself drive the slow structural re- 
laxation process at the time scale of the imposed strain 
rate. For oscillatory measurements, the natural scale of 
strain rate is its amplitude 70 = 70 w. Thus, if a struc- 
tural relaxation process is present, its time scale t(to) 
should depend on the applied strain rate amplitude, 
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However, in traditional strain or frequency dependent 
rheological measurements, the amplitude of the strain 
rate varies for each data point for both frequency and 
strain dependent measurements. This makes it difficult 
to discern the effects of strain rate on the time scale r(jo) 
of a structural relaxation process. 

We circumvent this problem by performing frequency 
dependent rheological measurements at constant strain 
rate amplitude 70. This 'constant rate frequency sweep' 
is achieved by applying a strain amplitude that is in- 
versely proportional to the oscillation frequency, as 
shown in the inset of Fig. 2. We perform a series of such 
measurements with 70 ranging from 0.01 s" 1 to 5s _1 , 
shown in the main graph of Fig. 2. For comparison, the 
linear viscoelastic measurement (see Fig. 1) performed at 
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FIG. 2: (Color online) Frequency dependent storage G'(u>) 
(solid symbols) and loss modulus G"(uj) (open symbols), mea- 
sured in constant rate oscillatory measurements at different 
strain rate amplitudes 70 = 701^; 70 = 0.01 s _1 (grey circles), 
0.05 s^ 1 ( green diamonds), 0.5 s 1 (red triangles), 5s 1 (blue 
squares). The linear oscillatory measurement of G' and G" 
are shown by the black circles. Inset: Strain amplitude 70 
versus u corresponding to the constant rate measurements in 
the main graph. Top: Corresponding ratios of the third to 
the first harmonic response I(3u)/I(ui) as a function of u>. 



a constant strain amplitude (circles) is also included. 

The viscoelastic response depends strongly on the 
strain rate amplitude. With increasing 70, the entire fre- 
quency dependent response is shifted towards higher uj, 
as expected from Eq. (1). Remarkably, however, the gen- 
eral shape of the frequency dependent storage and loss 
moduli is surprisingly insensitive to 70. At low frequen- 
cies the response is liquid-like: G" is larger than G' and 
both moduli scale with frequency as a power law, with 
a ratio of sa 2 between the corresponding exponents. At 
higher frequencies there is a pronounced peak in the loss 
modulus, followed by a shallow minimum and a final slow 
increase at the highest frequencies. The storage modu- 
lus rises continuously until it reaches a plateau where 
it shows only a very weak frequency dependence. The 
response remains surprisingly harmonic for all our mea- 
surements @ with a ratio of the third to the first har- 
monic component of the response that is always smaller 
than 20%; in the relevant frequency range around the 
peak position this ratio is smaller than 10%, as shown in 
the top graph of Fig. 2. This justifies our treatment of 
the data in terms of the viscoelastic moduli G' and G". 

To highlight the striking similarities in the shape 
of the viscoelastic response, we rescale the moduli in 
both magnitude and frequency, using G sca i e d*(w) = 
G* (uj /b(j )) /a(jo) , where 0(70) and 6(70) are the seal- 
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FIG. 3: Constant rate frequency sweep measurements shifted 
onto a single master curve. Inset: Corresponding amplitude 
and frequency shift factors 0(70) and 6(70) versus 70. 



ing factors for the magnitude and frequency, respectively. 
All the data can be superimposed onto a single master 
curve, as shown in Fig. 3. The scaling factors are plot- 
ted in the inset to Fig. 3 as a function of 70. The value 
of a(7o) varies by less than a factor of 3 over the entire 
range of 70 accessed, whereas the frequency scaling fac- 
tor 6(70) depends strongly on 70. At large strain rates 
we observe a power-law dependence, 6(70) oc 7q, with 
v w 0.9; at low strain rates it asymptotically approaches 
unity, which suggests that in this regime the imposed 
strain rate is sufficiently low that it no longer impacts 
the viscoelastic behavior. The rate dependence of the 
scaling factors is consistent with the behavior expected 
from Eq. (1), with 6(70) = tq/t^q). 

The observed scaling behavior is a direct consequence 
of a structural relaxation process, which leads to a de- 
crease in G'(u>) and a concomitant peak in G"(cu) in lin- 
ear viscoelastic measurements. However, this structural 
relaxation occurs at too low a frequency to be directly 
observed in the data in Fig. 1(a); nevertheless the incip- 
ient behavior is observed in the linear viscoelasticity at 
low frequencies and is responsible for the peak observed 
in G" in the nonlinear behavior. 

The structural relaxation time decreases as 70 in- 
creases, as expected from Eq. (1). The relaxation peaks 
in the data shown in Fig. 3 occur when 2tt/lo becomes 
comparable to t(t"o). To isolate this contribution due to 
structural relaxation, we must separate the peak from the 
subsequent rise in the data at higher frequencies. Closer 



FIG. 4: (a) Loss modulus G" for a series of constant strain 
rate oscillatory measurements with a high frequency contribu- 
tion Cy/uj subtracted (dotted line), (b) The uppermost curve 
shows the same data shifted onto a single master curve. Ad- 
ditional master curves shown: 80vol% oil in water emulsion, 
aqueous foam (Gillette Foamy), PMMA hard sphere suspen- 
sion at <f> m 0.58. Each data set has been scaled for clarity. 



inspection of the data reveals systematic differences be- 
tween the scaled curves at the highest frequencies. This 
suggests that a significant contribution to the response 
does not depend on strain rate at sufficiently high lu. 
The high frequency viscoelastic response in compressed 
emulsions exhibits a y^-dependence, attributed to vis- 
cous flow along randomly oriented slip planes (l6| : this 
results in a contribution that scales as y/tJ. The high 
frequency linear response of our system shows the same 
behavior, as shown by the dashed line in Fig. 1(a). We 
thus subtract a component c^/lu from each of the data 
sets, where the constant c is determined from a fit to the 
high lu response; a single value of c is used for all the data, 
independent of 7. The resultant loss moduli now exhibit 
well defined peaks, as shown in Fig. 4(a); the component 
subtracted is also shown as the dashed line. 

The shapes of the resulting relaxation peaks are re- 
markably similar at each value of 70. To highlight this, 
we scale all the data onto a single master curve by rescal- 
ing the frequency, as shown by the uppermost data set 
in Fig. 4(b); because we have subtracted the y/ui compo- 
nent, we do not need to scale the amplitude. The scaling 
of the data suggests that the shape of the relaxation spec- 
trum remains unchanged as the strain rate amplitude is 
varied. This should provide important new insight about 
the nature of the relaxation and the distribution of time 
scales that make up the peak. 

The observed scaling suggests that G* can be expressed 
as a function of frequency and strain amplitude: 

G>, 70) « G* R (cu/b( lo cu)) + G* HF (lu) , (2) 

where Gj^ (u>) describes the functional form of the lin- 
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ear viscoelastic behavior due to the structural relaxation 
process, and G^ F (u) describes high frequency contri- 
butions that are independent of 70; caution must be 
used for u < w max , where the slowest relaxation must 



dominate [17]. This expression for G* reflects both the 
linear and the incipient nonlinear response of the system. 
For low enough values of strain, Eq. (2) describes the 
complex modulus in the linear viscoelastic regime, with 
6(70) ~ 1. At larger strains, the structural relaxation 
shifts towards higher frequencies, scaling as 6(70). This 
allows us to establish a direct link between linear and 
nonlinear rheology through the strain-rate dependence 
of the structural relaxation time. To further illustrate 
the physical origin of this link, we provide a schematic 
movie in the supplemental information [3]. 

Our results provide an important new measurement 
technique for probing the rheological properties of soft 
materials. By maintaining a constant strain-rate ampli- 
tude during our measurement, we can isolate the struc- 
tural relaxation; by varying the amplitude of the strain 
rate, we force the frequency of this structural relaxation 
into a regime accessible to the rheometer. This allows 
us to directly probe this structural relaxation, even if 
it is at a frequency too low to be easily accessible. By 
scaling the data, we highlight the low-frequency contri- 
butions of the structural relaxation. A wide variety of 
soft materials exhibit such a low-frequency structural re- 
laxation; thus, this method should be very broadly ap- 
plicable. To illustrate this generality, we perform simi- 
lar experiments on several representative samples of soft 
materials, and isolate the structural relaxation by scal- 
ing the peak data onto master curves. Similar behavior is 
observed for all samples as shown in Fig. 4(b), where we 
plot master curves obtained for a compressed emulsion 
(80% silicon oil drops of w 50 /im dia in water, stabi- 
lized by 10 mM SDS), a foam (Gillette Foamy), and a 
hard-sphere colloidal suspension (PMMA in cycloheptyl- 
bromide/decalin). In all cases, good scaling is observed; 
moreover, differences in shape are also apparent, reflect- 
ing different contributions to the relaxation spectrum. 

The scaling employed here is reminiscent of an ap- 
proach that has been successfully employed for polymer 
melts. Rheological measurements of G'(oj) and G"(u>) 
for different temperatures are shifted onto a single mas- 
ter curve that reflects the viscoelastic behavior in a dra- 
matically extended range of frequencies. This technique 
is called time-temperature superposition (TTS) lij 2oj ]. 
By analogy, we call our approach strain-rate frequency 
superposition (SRFS). With SRFS, we can obtain de- 
tailed information on the strain rate dependence of the 
slow relaxation process. Our picture combines linear 
and nonlinear oscillatory rheology into a unified picture, 
thereby accounting for the unusual behavior typically ob- 
served in many soft materials. Our experimental results 
show that nonlinear viscoelastic measurements contain 
useful information on the slow relaxation dynamics of 



these systems. Their yielding behavior directly probes 
the structural relaxation process itself, shifted towards 
higher frequencies by an applied strain rate. Physically, 
this suggests that this structural relaxation, driven by 
an imposed strain rate, results in the same response as 
the equilibrium structural relaxation at a much lower fre- 
quency. This is a surprising result, as for the steady shear 
case theory[T3] suggests a rate-dependent change of the 
high frequency flank of the relaxation; for the systems 
studied we do not observe this in our oscillatory mea- 
surements. By exploiting this behavior, SRFS isolates 
the shape of the structural relaxation, even if it is too 
slow to be accessed with conventional oscillatory rheol- 
ogy. Thus, SRFS will provide new insight into the phys- 
ical mechanisms that govern the viscoelastic response of 
a wide range of soft materials. 
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